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A new formulation for coUisional kinetic theory is presented based on the use of Lie- 
transform methods to eliminate fast orbital time scales from a general bilinear collision 
operator. As an application of this new formalism, a general guiding-center bilinear Fokker- 
Planck (FP) collision operator is derived following the elimination of the fast gyromotion 
time scale of a charged particle moving in a nonuniform magnetic field. It is expected 
that classical transport processes in a strongly magnetized nonuniform plasma can, thus, 
be described in terms of this reduced guiding-center FP kinetic theory. The present paper 
introduces the reduced-collision formalism only, while its applications are left to future 
work. 

PACS numbers: 52.25.Dg, 52.25.Fi, 52.20.Dq 



1 



I. INTRODUCTION 



Lie-transform perturbation methods are normally associated with the asymptotic elim- 
ination of fast time scales from single-particle Hamilton equations d 12] • In the present 
paper, the general Lie-transform rules for the transformation of an arbitrary bilinear colli- 
sion operator are presented. 

As a first step, Lie-transform methods are used to asymptotically eliminate a fast orbital 
time scale from the coUisionless (Vlasov) evolution operator of a coUisional kinetic equation. 
This dynamical reduction is performed through a phase-space coordinate transformation, 
from which coUisionless reduced Vlasov kinetic theories are derived (e.g., gyrokinetic Vlasov 
theory As a result of this coordinate transformation, the collision operator is trans- 

formed into a new collision operator in which the fast time-scale dependence is still present; 
this feature is common to all coUisional theories based on the use reduced guiding-center 
coordinates jSUHllZl- To eliminate this residual fast-time-scale dependence, we proceed with 
a second step involving a fast-angle averaging procedure and a closure scheme based on a 
time-scale ordering involving the characteristic coUisional time scale. A reduced coUisional 
kinetic equation is, therefore, obtained in which the fast-angle-averaged distribution evolves 
in a reduced phase space. Note that fast-angle averaging is required only in the presence 
of collisions since the first step involves a reversible phase-space coordinate transformation 
that requires no averaging. 

The present approach is to be contrasted with the traditional multiple time-scale ap- 
proach [H], whereby the coUisional kinetic equation in particle phase space is expanded 
in powers of a small parameter e, with the particle distribution function expanded as 
/ = I]n=o fn- In general, the zeroth-order component /o is independent of the fast 
time scale and the fast-time-scale dependence of fn involves contributions from coUision- 
less (Hamiltonian) dynamics and coUisional dynamics, which introduces several layers of 
tedious algebra before a reduced coUisional kinetic equation is derived [3 CHI, with a lin- 
earized collision operator. 

By using Lie-transform methods, we obtain simpler and more compact expressions 
for transformed guiding-center Fokker-Planck (FP) collision operators (for uniform and 
nonuniform magnetic fields) when compared to those obtained by the standard approach 
for uniform magnetic fields Ej. The present approach is expected to be appropriate 
appropriate for applications in coUisional gyrokinetic theory and gyrokinetic particle sim- 
ulations [THIT^. 

A. Brief Introduction to CoUisional Plasma Kinetic Theory 

We begin with a brief introduction of coUisional plasma kinetic theory (see Refs. |13j - 
|15j . for example). The study of irreversible transport processes in nonuniform plasmas is 
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based on solutions of a collisional plasma kinetic equation 



dt 



(z,t) = C[/](z,t) 



(1) 



where / is the test-particle distribution function and z are coordinates in the six-dimensional 
single-particle phase space. The Vlasov operator 



characterizes the dissipationless time evolution of /, where the Hamiltonian particle orbit 
z(t; zq) is the solution of the Hamilton equations z'^{z,t) = {z", h{z,t)} with the initial 
condition zq; here, h and { , } are the single-particle Hamiltonian and Poisson bracket, 
respectively. (Summation over repeated indices is implied throughout the paper; greek 
indices take values from 1 to 6, while latin indices take values from 1 to 3.) 

In the absence of collisions (C = 0), Eq. ((T)) is the Vlasov equation [12]: df{z, t)/dt = 0, 
whose solution /(z, t) is constant along a Hamiltonian particle orbit. The collision operator 



in Eq. (0) characterizes the dissipative (irreversible) time evolution of /, where the bi- 
linear collision operator C[/; /'] describes binary collisions between test-particles and field- 
particles (with distribution /') ; the summation is over all field-particle species (and includes 
like-particle collisions) . 

Solutions of the collisional kinetic equation (P) are quite difficult to obtain in general 
and various approximation schemes must be adopted in order to arrive at useful solutions. 
One such approximation scheme involves removing fast orbital time scales associated with 
Hamiltonian particle orbits in six-dimensional phase space E|- Here, the motion of 
magnetically-confined charged particles exhibits three distinct orbital time scales (r^ -C 
Tb -C Td): (1) a fast gyration time scale Tg associated with the gyromotion of charged 
particles about a magnetic field line; (2) an intermediate bounce time scale Tb associated 
with the periodic parallel motion of charged particle along a magnetic field line; and (3) 
a slow drift time scale associated with the perpendicular motion of charged particles 
across nonuniform magnetic field lines. 

After choosing a slow time scale of interest (labeled r) , we identify all fast time scales 
which satisfy the condition -C r. For example, if we are interested in plasma dynamics on 
the bounce time scale (r ~ u) then the gyromotion time scale is considered fast (r^ <^ r) 
and it can be asymptotically removed. For each such fast orbital time scale r^, a pair 
of action-angle variables ( J^, 6'') is assigned: 6'' is the fast angle variable (with 6'' = 
denoting the fast orbital frequency) and its canonically-conjugate action Jk is an adiabatic 
invariant on a time scale longer than = § dd^ /VL^ (i.e., the fast-angle average of Jk is 
zero). The asymptotic elimination of a fast orbital time scale from the Vlasov operator 




(2) 



(3) 



3 



(2) can be carried out to arbitrary order in the small parameter = uo/Vt^ ^ 1 (where 
idr oj denotes a characteristic frequency) jH HHj . 

In the present paper, we focus our attention on the asymptotic elimination of the gy- 
romotion time scale leading to the reduced guiding-center dynamics of charged particles in 
nonuniform magnetic fields. Here, the reduced guiding-center coordinates (K,S, fi,6) in- 
clude the guiding-center position X, the guiding-center kinetic energy £, the guiding-center 
magnetic moment n, and the guiding-center gyro-angle 6. Each coordinate is expressed as 
an asymptotic expansion in powers of the dimensionless parameter = p/Lb 11211201121], 
where p is the characteristic gyroradius of a charged particle and Lb denotes the magnetic 
field nonuniformity length scale. 

B. Lie-transform Phase-space Transformations 

The standard mathetical method used to asymptotically eliminate fast time scales is 
the Lie-transform method, which we briefly review here. In particular, we show the trans- 
formation rules by which each element of the collisional kinetic equation (P) (i.e., /, d/dt, 
and C) transform under a near-identity transformation. 

Using Lie-transform methods PP, the process by which a fast time scale is removed from 
Eq. Q involves a near-identity transformation on particle phase space: 

%: z ^ Z(z; e) = %z, with Z(z; 0) = z, (4) 

where e -C 1 is a small parameter (for reduced guiding-center dynamics, e = cb)- In what 
follows, the new six-dimensional phase-space coordinates Z = (Zr; J, 6) are divided into 
the reduced phase-space coordinates Zr (the reduced particle dynamics takes place in a 
/our-dimensional space) and the fast canonical action-angle pair {J, 6). 

Next, using the transformation (jH), we define the push-forward opeiatoi on scalar fields 
P] induced hj the near-identity transformation Pj): 

T;': f ^ F ^ T7V, (5) 

i.e., transforms a scalar field / on the phase space with coordinates z into a scalar field 
F on the phase space with coordinates Z: -F(Z) = T7^/(Z) = f{T~^Z) = /(z). Since the 
transformation is invertible, i.e., there exists an inverse near-identity transformation 

T^'^ : Z ^ z(Z;e) = T'^Z, with z(Z; 0) = Z, (6) 

we also define the pull-back opeiatoi jT]: 

T, : F ^ / = T,F, (7) 

i.e., transforms a scalar field F on the phase space with coordinates Z into a scalar field 
/ on the phase space with coordinates z: /(z) = T^F{z) = F{%z) = F{Z). 
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Lastly, we consider the transformation of a general operator C : f ^ Cf acting locally 
on the particle distribution / such that Cf is also a scalar field on particle phase space. 
Hence, the induced transformation of the scalar field Cf follows the push-forward rule 
(0): Cf — > T~^Cf. By combining this push-forward transformation with the pull-back 
transformation / = T^F, we construct a new local operator C^ : F C^F defined as 



This induced transformation rule for operators will be used in the next Section to transform 
the Vlasov operator d/dt and the collision operator C appearing in the collisional kinetic 
equation ((T)). 

C. Organization 

The remainder of the paper is organized as follows. In Sec. II, we investigate the fast- 
time dependence of the transformed collisional kinetic equation and discover that, although 
the fast time scale has been eliminated from the dissipationless (Vlasov) part of the kinetic 
equation, the dissipative (collisional) part has retained fast-time dependence as a result of 
collisions. This is not too surprising after all since these phase-space transformations were 
never designed for that purpose. This means that the transformed collision operator is not 
yet reduced in the sense that the fast-angle-independent (averaged) part, denoted (F), and 
the fast-angle-dependent part, denoted F = F — (F), of the transformed distribution func- 
tion F are still coupled by collisions. The reduction of the transformed collision operator 
involves a closure scheme yielding a reduced collision operator acting on (F) alone. This 
closure scheme is based on an asymptotic expansion based on a new expansion parameter 
= u /VL = pI Xi, <^ 1 defined as the ratio of the characteristic collision frequency z/ over 
the fast frequency Vt or the ratio of the gyroradius p over the collisional mean-free-path 
For most practical applications, however, we argue that only the lowest term in this 
ejy-expansion is relevant. 

The purpose of the present paper is, therefore, to derive first-order (e^ = p/Lb) cor- 
rections of a reduced guiding-center Fokker-Planck (FP) collision operator, which might be 
suitable for analytical and numerical study of collisional and turbulent transport processes 
in nonuniform magnetized plasmas. In Sec. Ill, we derive the lowest-order term in a reduced 
Fokker-Planck (FP) collision operator. Starting from a standard FP collision operator (de- 
scribing collisional drag and diffusion in momentum space), we derive an expression for a 
reduced FP collision operator in a reduced (4 + 1) phase space with coordinates (Zj:j, J). In 
particular, the reduced FP collision operator describes spatial drag and diffusion processes 
normally described from a collisional fiuid-moment approach. 

As a specific example, we focus our attention in Sec. IV on the guiding-center phase- 
space transformation associated with the asymptotic elimination of the fast gyro-motion 
time scale of a charged particle moving in a nonuniform magnetic field. General expressions 
are obtained for the collisional FP drag and diffusion coefficients in guiding-center phase 




(8) 



5 



space. For clarity of exposition, we also use an isotropic field-particle model introduced 
in Sec. II to investigate how magnetic-field nonuniformity influences these guiding-center 
FP coefficients. We also discuss possible extensions of our work introduced by background 
nonuniform electric fields (which are outside the scope of the present work) and summarize 
our findings in Sec. V. Lastly, we present some technical details about the guiding- center 
transformation for nonuniform magnetic fields in Appendix A. 



II. TRANSFORMED COLLISIONAL KINETIC EQUATION 



The transformation of the coUisional kinetic equation (^J), induced by the near-identity 
phase-space transformation Q, yields a new collisional kinetic equation: 

|F(Z,t) = CjF](Z,t) = i: C^,,^[F;F']{Z,t), (9) 

where F = T^^f is the transformed test-particle distribution function obtained from the 
push-forward of the test-particle distribution function /, the transformed Vlasov operator 
d^/dt is defined from Eq. (jH)) as 



|^,z,,).t-(|t,f)(z,o, 



(10) 



while the transformed bilinear collision operator is defined as 

C(,,,,)[F;F'](Z,t) ^ T:\C[T,F;T,F']){Z,t), (11) 

where e' = u/Q' is the small parameter associated with the field-particle species. We 
note that the bilinearity of the collision operator Q allows a separate treatment of the 
test-particle (/) and field-particle (/') distributions and, for convenience of notation, we, 
henceforth, write the transformed collision operator (lllj) as 

C,[F](Z,t) ^ T7i(C[T,F])(Z,t), (12) 

and omit mention of the field-particle species unless otherwise needed. 

The transformed Vlasov operator ^TU^ is now explicitly written in terms of the new 
phase-space coordinates (Z/?; J, 9) as 

_ -9 ■ 9 _ dK ■ d 
Jt - dl^ ^^dZ^ - Tt ^ ^^^^ 

where the transformed Hamilton equations Zf = {Z°', H}^ are defined in terms of the 
transformed Hamiltonian H = T^^h and the transformed Poisson bracket 

{F, GUZ) ^ T:'{{T.F, T.G})(Z). (14) 
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We note that the transformed Hamiltonian H{Zji; J) and the transformed Poisson bracket 
{ , }e are (by construction) independent of the fast angle 9, so that Hamilton's equations 
Zr(Zr; J) and 9 = ^(Zr; J) - and, thus, the transformed Vlasov operator d^/dt - are 
independent of the fast angle 9 while J = (to all orders in e). 

The 6'-independence of the transformed Vlasov operator (fT^ allows the separation of 
the transformed collisional kinetic equation Q into two coupled kinetic equations: one for 
the ^-avera^ed distribution, denoted (F), and one for the 6'-de^>en(ient distribution, denoted 
F = F — (F). The collisional kinetic equation for (F) is obtained by ^-averaging both sides 
of Eq. ©: 

^(F) = (C,[F]) = (C,[(F)]) + (^[F]), (15) 

where we used the 6'-independence of the reduced Vlasov operator d^/dt on the left side 
of Eq. (jl5|) . while the collisional kinetic equation for F is obtained by subtracting Eq. (jl5|) 
from Eq. Q: 

Note that the term {Cf\F]) in Eq. ()15|) does not vanish in general and, thus, the collisional 
kinetic equations (fT3j) and (fTBj) are coupled through collisions (i.e., when C = 0, it can be 
shown that F = to all orders in e). Hence, Eq. (fT3j) still contains a dependence on the 
fast angle 9 through the ^-dependent part F of the transformed distribution F. 

The ultimate goal in reduced collisional kinetic theory is to transform Eq. |TH|l into a 
closed collisional kinetic equation: 

^(F)(ZR,t;J) ^ CR[(F)](ZR,t;J), (17) 

where the reduced (R) collision operator Cr acts on the ^-independent part (F) alone. To 
obtain this collision operator, we must solve the collisional kinetic equation (fTH|) for F in 
terms of (F) and substitute this functional solution into (Ce[F]). This is accomplished by 
expanding the transformed distribution function F as a Fourier series in the fast orbital 
angle 9: 

F = {F) + [Fee''' + F/ e"' 
e>o 

where F/ = F_^, and defining the averaged collision operators 



dig] = {e-'''C,[ge'''\), (19) 

where g is an arbitrary ^-independent function and {k, i) are arbitrary integers. By sub- 
stituting these expressions into Eqs. (fT^ - (fTH|) . respectively, we obtain 

^{F) = C°[(F)] + EC^[^^], (20) 

(ziQ + ^]F, = C°[(F)] + EQl^^]- (21) 
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The approximate solution of Eq. (|2ip for each integer £ 7^ is based on the small parameter 
= V /Vt <^ 1 (which is different from e = uj/VL), where v is the characteristic collisional 
dissipation rate. To first order in e^, and assuming that C°[(F)] 7^ in Eq. (PTj) . we, 
thus, find the functional solution F(,[{F)] = C^[{F)], which when substituted into 

the second term on the right side of Eq. (|^n|) yields the first two terms in the asymptotic 
expansion for Cr: 



Cr[{F)] = C°[(F)] - ^e.J:C^\^ ] + O (ee„ el) , (22) 

where 



C°[(F)] ^ (C,[(F)]) = {T:'C[TAF)]) (23) 

is the lowest-order term in the closed transformed collision operator. 

The reduced collisional kinetic equation (fTTjl . including the reduced collision operator 
()22|l . is obtained from the collisional kinetic equation as a result of two asymptotic 
expansions in powers of e and e^- In most practical applications, however, the zeroth-order 
expression is sufficient and, thus, we, henceforth, focus on applications of this lowest- 
order term in the development of reduced collisional kinetic theories. We note that the 
two small parameters and introduced above serve very different purposes: = p/Lb 
is involved in the elimination of fast gyro-motion time scale from the Vlasov operator 
and = u/Q = p/X^ is involved in the elimination of fast gyro- motion time scale from 
the collision operator. The ratio of these two dimensionless parameters introduces a new 
parameter A = es/ty = X^/Lb in the collisional kinetic theory of irreversible transport 
processes in nonuniform magnetized plasmas jT3]. According to neoclassical collisional 
kinetic theory [Hj, for example, the c/asszca/ collisional regime is identified by the condition 
A <^ 1, i.e., the collisional mean-free-path is much shorter than the magnetic nonuniformity 
length scale, so that magnetic nonuniformity can be ignored at the lowest order in e^. In 
the neoclassical legmve, on the other hand, we have A > 1 (i.e., Xy > Lb), while in the 
collisionless (long-mean-free-path) regime, we have A ^ 1 (i.e., X^ ^ Lb)- In the latter 
two cases, the motion of a guiding-center particle explores a significant portion of the 
magnetic field nonuniformity before it encounters a collision. 



III. TRANSFORMED FOKKER-PLANCK COLLISION OPERATOR 



The collisional operator C has so far been treated as an arbitrary bilinear collision 
operator. We now focus our attention on an important class of bilinear collision operators 
used in plasma physics: the Fokker-Planck (FP) collision operators (see, for example, 
Refs. ^Hl 1221 ) • The general form of the FP collision operator is expressed as 



Cfp[/](x,p) 



d_ 



i^^(x,p)/(x,p) - D^^(x,p) 



dpi 



(24) 
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where z = (x, p) are phase-space coordinates for a test-particle of mass m and charge 
e (here, x denotes the particle position and p = mv its kinetic momentum); the time 
dependence is suppressed for simplicity of notation and the momentum representation is 
used to facilitate the analysis that follows. The collision operator describes momentum- 
space scattering events in both unmagnetized plasmas [2^ and magnetized plasmas [^ 1^: 
the three-component vector K = X]' describes coUisional momentum drag while the 

three-by-three symmetric matrix D = J2' D[/'] describes collisional momentum diffusion; 
the coefficients K[/'] and D[/'] are functionals of the field-particle distribution function /'. 

For small-angle-deflection Coulomb collisions in an unmagnetized plasma involving test- 
particle species a and field-particle species b [Hj, for example, the coefficients Kafe[/fc] and 
Dabi/fe] are expressed in terms of Rosenbluth potentials {Hab,Gab) as 

K„,(z) = f fr..^ /.VAM^^I^:^) , 

op \ mi) J |v — v'l y op 

(25) 

^^"^'^ = 2aF9Fr"'^"i^'^^^'^^^"""^'"""'J = 2^F^' 

where Tab = 47re^e^ In A and dp • Dab = i^b/^a) ^ab- Note that magnetized FP collision 
operators can also be expressed in terms of Rosenbluth potentials However, the 

unmagnetized FP collision operator ()24|) can still be used to study collisional transport 
processes in magnetized plasmas if the characteristic gyroradius is larger than the Debye 
length for each plasma particle species (i.e., the Alfven speed is sub-luminal for each plasma 
particle species), which is well satisfied for electrons and ions in high-temperature tokamak 
plasmas. 

When the field-particle distribution is isotropic in momentum space such that the 
Rosenbluth potentials -ffafe(z) = Habi'x, p) and Gabi'z) = Gab(x, p) depend only on the 
magnitude of the particle momentum, the FP coefficients become 

Ka6 = -VabP and Dab = Diiab ^ + D±ab (l - ^ j , (26) 

where the collisional momentum-drag frequency is I'ab = — i)- / v) H'abip) , the diffusion 
coefficients parallel and perpendicular to p, respectively, are D\\ab = (l/2)G^^(p) and 
D±ab = (l/2p) G'abip)- The isotropic field-particle model will become useful later on when 
we investigate the effects of magnetic-field nonuniformity on guiding-center collisions. 



A. Standard Transformation Procedure 



The standard transformation procedure for the FP collision operator 
in terms of the operator transformation Cpp — ^ Ctfp'- 



Ctfp[F'] 



1 A. 

J 



OF ' 



is expressed 



(27) 
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where J denotes the Jacobian of the general transformation z = (x, p) — Z = Tz, 
the distribution function F = T~^f denotes the push-forward of the particle distribution 
function / generated by the general transformation T, and the new FP coefficients in 
Eq. are defined as 

where we have used standard properties of the push-forward operator T~^: T~^(/ g) = 
(T~^/) {T~^g). Note that, because the transformed FP collision operator ()27p is written in 
divergence form, it is guaranteed to satisfy the particle conservation property. 

Although the standard transformation procedure (j27j) is straightforward to apply for fi- 
nite transformations y/^ (e.g., transformation from Cartesian to cylindrical momentum co- 
ordinates), the explicit calculation of the transformation-matrix components T~^{dZf / dp^) 
can be difficult to implement in practice for near-identity transformations on six-dimensional 
phase space, such as those contemplated in the present work. The difficulty resides in the 
fact that theterms T~^{dZ'^ / dp^) in Eqs. (f^Hj) - (pi?jl involve a double asymptotic expansion 
in powers of e El IZ| • To bypass this difficulty, we introduce a new formulation of the FP 
collision operator based on the use of the noncanonical Poisson bracket { , } used in 
describing single-particle Hamiltonian dynamics. 



B. Poisson-Bracket Formulation 



First, we point out that the general FP collision operator ()24|) can be written as 

Cfp[/](z) = - [x\ [k^(z)/(z) - D^^(z) {x^ /(z)}] }, (30) 
in terms of the noncanomca/ Poisson bracket 
9/(z) dg{^) 



(9p 



9/(z) dg{'z) ^ eB(x) 



9p 



(9x 



dp dp 



(31) 



where / and g are two arbitrary phase-space functions and the magnetic field B appears 
in Eq. ()31|1 because p = mv is not the canonical momentum. 

Next, using the transformed Poisson bracket ()14|). the transformed FP collision operator 
is obtained from Eq. (fT^ as 



C,pp[F](Z) = - {X:(Z), [kI{Z)F{Z) - D^iZ) {XKZ), F(Z)}J } 



(32) 



where the transformed FP coefficients are Ki{Z) = J-^K\Z) and Di^{Z) = T;^D'^{Z), 
and Xl{Z) = T~^x^{Z) = represents the particle position x expressed as a function of 
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the transformed phase-space coordinates x = Xe(Z), and X represents the transformed 
spatial coordinate. 

We now point out that a general Poisson bracket in six-dimensional phase space (with 
coordinates Z") can be expressed in terms of a six-by-six antisymmetric Poisson matrix 
Jf^ = {Z"', Z^je and the Jacobian JT" as a phase-space divergence 



dF dG 



1 d 



dZ 



dzp 



J dZ^ 



,^ dZP ^ J 



where the functions F and G are arbitrary functions on the transformed phase space, while 
the Jacobian J' and the anti-symmetric Poisson matrix J"^ satisfy the Liouville identities 
da{<J Je^) = (for (3 = 1, ...6). With this Poisson-bracket formulation, we may rewrite the 
Fokker-Planck collision operator ()32|) in the same phase-space divergence form as Eq. ()27|): 



C.FP[F](Z) 



1 d 



al3 



dF 



JdZ^ ^ ^ ' dZP 
where the phase-space coUisional drag vector is defined as 

while the phase-space collisional diffusion tensor is defined as 

Vf = - Jf 9^X, • D, • d^X, Jf = • D, 



(33) 



(34) 



(35) 



which is also symmetric since is symmetric. In Eqs. fl34j) and we have introduced 
the vector-valued components A" defined as 



al3 



dZ^'' 



(36) 



which offer computational advantages over the standard transformation approach (de- 
scribed above) since asymptotic expressions for the antisymmetric Poisson matrix ^ and 
Xg are readily available for the guiding-center \20\ and bounce-center |25| near-identity 
phase-space transformations. 

In the next Section, we derive expressions for the guiding-center Fokker-Planck collision 
operator based on the general phase-space divergence form dSBI), with coefficients given by 
Eqs. dSl and dnH). 



IV. REDUCED GUIDING-CENTER FOKKER-PLANCK OPERATOR 



As an application of the Poisson-bracket formalism presented in Sec. III.B, we consider 
the transformation of a general FP collision operator under the guiding-center transforma- 
tion pUJ: z = (x, p) ^ Z = (X,S] fi,9), where X denotes the guiding-center position, S 
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denotes the guiding-center kinetic energy {S = ]5|/2m + /i-B), /i denotes the guiding-center 
magnetic moment, and 6 denotes the guiding-center gyroangle. The choice of the actual 
velocity-space guiding-center coordinates used in the guiding-center FP collision operator 
can vary depending on its applications. Here, we use the {S,fi) coordinates because of the 
simplicity of the derived guiding-center FP collision operator. Future work will consider 
replacing the magnetic moment fi with the guiding-center's pitch angle = arccos(p||/p). 

A. Nonuniform Magnetic Field 

In the collisionless long-mean-free-path regime A = es/^u = ^u/Lb ^ 1 (where the 
collisional mean-free-path Aj, is much longer than the magnetic nonuniformity length scale 
Lb), the guiding-center motion of a charged particle samples a significant portion of the 
magnetic-field nonuniformity before it encounters a collision and, thus, magnetic field 
nonuniformity becomes important when describing small-angle collisions. 

According to the guiding-center transformation [2011^5 the particle position 

X = X + = X, (37) 

is expressed in terms of the guiding-center position X and the gyroradius vector (whose 
expression includes corrections due to magnetic- field nonuniformity); the dimensionless 
parameter e = es represents the magnetic field nonuniformity. The expression for the 
gyroradius vector (up to order e) 

P.(Z) = Po - e(Gf + ^Gi-rfpo) + ••• (38) 

is given in terms of the components of the first and second order generating vector fields 
(G"j*, G2) for the guiding-center transformation [2111 E] and contains the lowest-order gyro- 
radius vector Pq = — Gf- as well as the first-order correction p^ = — — \ Gi ■ d p^, an 
explicit expression for the gyroradius first-order correction p^ is not needed in this Section 
and is given in Appendix A [see Eq. ()A.5|) ]. 

The Jacobian J = mB^^/\v\\\ for the guiding-center transformation [20] is defined in 
terms of the guiding-center phase-space function -B|| (X, S, fi) = h-'B* derived from the gen- 
eralized magnetic field B* = B + e {cp\\/e) V X b, where |f||| = \p\\\/m = ^J (2/m) {£ — fiB). 
Next, the guiding-center Poisson bracket of two arbitrary functions F and G of (X, S, fi,6) 
is 

_n\dF(dG dG\ (OF dF\dG- 

V*i^ ^ - ^ V*Gj - ^ • V*F X V*G, (39) 
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where V* is a gradient operator defined in Appendix A [see Eq. ()A.4|) ] and the gyro-averaged 
guiding-center velocity is 



Vgc = v\\h + — ^ X \ fiVB + mvf b • Vb ) , (40) 

with n| = e-By/(mc) = Q + v\\ b • V X b. Lastly, the particle momentum Pe(X, 
T^^p is written in guiding-center coordinates as 



m — ; — = m 
dt 



Vgc + ( ^ ^ + Vge • V* 1 



(41) 



An explicit expression for is given in Appendix A [see Eq. ()A.6|) ]: here, we simply note 
that Pe satisfies the following properties 

I Pep = IpP = 2m £ and (p^) = mvgc, (42) 

up to first order in e. 

B. Guiding-center Fokker-Planck Collision Operator 

We now proceed with obtaining general expressions for the phase-space collisional drag 
vector and the phase-space collisional diffusion tensor associated with the guiding- 
center transformation (x, p) (X.,S, fi,9). First, using the definition (jHBj) for A"^ and 
the guiding-center Poisson bracket ()39|) . the guiding-center vector- valued components of 
A^, = {X„ are 

Age = -"-^b.dlX, ^^^^ 
< = + {p., ^}ge = |, (44) 



A'^ = (45) 



and we have ignored the component-vector Ag^ since it is not needed in evaluating the 
lowest-order guiding-center Fokker-Planck collision operator (j23j) . 

With these component-vectors (a greek-letter index, henceforth, excludes the gyroangle 
component), the guiding-center phase-space collisional drag vector and the guiding-center 
phase-space collisional diffusion tensor have the following guiding-center phase-space com- 
ponents 

/Cg", = K,. {X„ Z^} and Pgf = - {Z", X,} ■ D, • {x„ Z^} . (46) 

V. J gc 
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The lowest-order guiding-center Fokker-Planck collision operator is, thus, expressed in 
terms of the gyroangle-averaged guiding-center distribution (F) as 



J ( {K^^) (F) - (Dgf ) 



^\-^ / 



(47) 



gc / 

where expressions for the guiding-center FP coefficients {K^^) and {D^j?) are given with 



terms up to order e^. Note that the guiding-center Fokker-Planck collision operator now 
contains spatial derivatives and, therefore, Eq. ()47j] describes spatial coUisional drag and 
spatial coUisional diffusion. 

1. Guiding- center coUisional drag vector 

The gyroangle-averaged spatial components of the guiding-center coUisional drag vector 
needed for the lowest-order guiding-center FP collision operator (j23|) are 

where the omitted terms are of second order in e^. The gyroangle-averaged energy and 
magnetic-moment components, on the other hand, are 

We note that the components P^ -ffSU j) of the guiding-center phase-space coUisional drag 
vector are non-vanishing in the absence of magnetic nonuniformity and, thus, magnetic 
nonuniformity adds only small first-order corrections to these components. We expect the 
spatial components (jlH|) of the guiding-center phase-space coUisional drag vector, on the 
other hand, to vanish in the absence of magnetic nonuniformity as will be demonstrated 
when we apply the isotropic- field-particle model to Eq. (jlSj) . 

2. Guiding- center coUisional diffusion tensor 

The gyroangle-averaged spatial components of the guiding-center diffusion tensor are 

where higher-order terms are omitted since spatial diffusion involves second-order spatial 
gradients already. The gyroangle-averaged energy and magnetic-moment diffusion coeffi- 
cients are 

Dl!) = (^-D.-Pi), (52) 
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- (St-.- St) • 

and the gyroangle-averaged mixed-component spatial diffusion coefficients are 

Once again we note tliat tlie components (f5T| - (f54|l of tlie guiding-center pliase-space col- 
lisional diffusion tensor are non-vanisliing in tlie absence of magnetic nonuniformity and, 
tlius, magnetic nonuniformity only adds small corrections to these components. We expect 
the spatial components (fS^ - lj^Uj) of the guiding-center phase-space collisional diffusion ten- 
sor, on the other hand, to vanish in the absence of magnetic nonuniformity. 



C. Isotropic Field-Particle Model 



To obtain a specific expression for a guiding-center FP collision operator (jTfj) and 
investigate how a nonuniform magnetic field affects the form of the FP collision operator, 
we now consider the case of an isotropic field-particle (IFP) distribution [see Eq. (pUj) ]. for 
which the transformed collisional drag vector and the collisional diffusion tensor are 



- l^e Pe, 



Dii 



2m£'' 



(57) 
(58) 



where z/^, Dy^, and D^^^ are calculated from Rosenbluth potentials and are assumed to 
depend only on the guiding-center energy £ (to lowest order in spatial nonuniformity) for 

we obtain the IFP identities 

(D||, p,) = D\\^2m8, (59) 



simplicity of presentation. From Eqs. (j37|) - (j3H|) and 

Ke • Pe = — 2m£ and p^ • • = p, 

and the gyroangle-averaged expressions 

(Ke) 
(De-Pe) 

b X (D,) X b 



gc-i 



Vs mv. 



(60) 



where 1 1 = I — b b and we introduced the convenient notation 



Dii 



m 



and V,, 



/i 



D\\e - Da 
2m S 



(61) 
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We now substitute Eqs. (f^ - (j3H|) and (jHUjl - ljUTI) to obtain explicit expressions for the 
guiding-center FP coefficients (pS |) - (j3m) . 

First, using the gyroangle-averaged expression for D^, the guiding-center FP spatial 
diffusion tensor becomes (to lowest order in e^) 

which is identical to the guiding-center spatial diffusion coefficient derived by Xu and 
Rosenbluth [TT]. Using the identities (|^. we also easily find, from Eqs. and the 
following guiding-center FP coefficients 



2£, (63) 



which are unaffected by magnetic-field nonuniformity as a result of the fact that the guiding- 
center and particle kinetic energies are equal to each other up to ffist order in [see 
Eq. (??)]. 

Next, using the gyroangle-averaged expressions (jEHI), the guiding-center FP spatial drag 
coefficient ()48|) and the guiding-center FP off-diagonal energy coefficient ()55|) become 

''^^ ^ f ^' U X ... (64) 



gc 

For the remaining guiding-center FP coefficients [ (j^Uj) . (jMj) . and we find 



(2 - es Age) ( Jf' /i, (65) 



D^!^) = (1 - Age) 2/i (2V, + ^ ) , (66) 
where Age = {v\\/Q) b • V X b denotes the guiding-center vorticity parameter, and 

D^'') = e^P^^^XYge, (67) 

where gyrorangle averages in Eqs. (jM|l - (jH7jl were computed using formulas found in Ap- 
pendix A [see Eqs. ()A.7j) - ()A.10j) ]. As was indicated above, the guiding-center FP spatial 
drag and diffusion coefficients (jMj) and ()67|) vanish in the limit of a uniform magnetic field 
(i.e., es = 0), while the guiding-center FP velocity-space drag and diffusion coefficients 
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(j65|) - (j66|) exhibit first-order corrections due to magnetic- field nonuniformity. 
D. Guiding-center Ordering and Beyond 

1. Guiding- center ordering 

In the standard guiding-center ordering, where <C 1, the first-order corrections in the 
IFP Fokker-Planck coefficients are indeed quite small and can, thus, be neglected. 

The non-vanishing contributions (jMj) and (jU7j) associated with spatial drag and diffusion, 
on the other hand, yield terms of order comparable to the IFP spatial diffusion term if 
magnetic and plasma nonuniformity length scales are comparable (i.e., p|V_Lln(F)| ~ e^) 
as might be expected in compact high-temperature tokamak plasmas. 

2. Background electric field 

We now discuss how background electric fields might be included in the IFP model flH^ - 
(jnZj), with separate orderings eg and e'^ associated with the background E x B velocity 
and its derivatives j^H 12^] , respectively. 

First, the term b X Vgc in the IFP spatial components (jMj) and ()67|) is replaced with 
b X (e^ Vgc + e^; ve), where eE denotes the ratio of the E x B velocity to the characteristic 
thermal speed (which is small in the core region of tokamak plasmas pi] but is larger in the 
edge region of tokamak plasmas with sheared flows [2^1 )• Next, the term Age in the velocity- 
space components is replaced with Age + (b/fi) • V X ve, where e'^ = p/Le 
denotes the ratio of the characteristic gyro-radius to the sheared-fiow length scale. Here, 
the second term describes the orbit-squeezing effect of edge electric fields |27j, which might 
lead to important modifications of collisional transport theory [221 12Z1 12H] ■ Future work 
will investigate the role played by background electric fields, which are relevant to the 
numerical gyrokinetic simulations of edge plasmas pUj . 

V. SUMMARY 

In summary, we have derived a general reduced Fokker-Planck kinetic equation by Lie- 
transform methods. The collisional kinetic equation (jl7j) was obtained as a result of two 
independent asymptotic expansions: the first expansion (with small parameter e = u/Q) 
was associated with the elimination of the fast time scale from the Vlasov operator; while 
the second expansion (with small parameter = u/Q) was associated with the elimination 
of the fast time scale from the transformed collision operator. 

As an application of the dynamical reduction formalism for bilinear collision operators, 
we derived a general guiding-center Fokker-Planck collision operator (j47p . with coefficients 
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given in compact form by Eqs. (j48|) - (j56|) . In order to show explicitly how magnetic-field 
nonuniformity affects the form of the guiding-center FP collision operator, we also made 
use of the simple isotropic-field-particle model (j^^ and obtained explicit expressions for the 
guiding-center FP coefficients ()62|) -f IU7j) . In this simplified model, the effects of magnetic- 
field nonuniformity are represented by the gyroangle-averaged guiding-center velocity Vgc 
and the guiding-center vorticity parameter Age- 

Future work will consider applications of the formalism leading to the reduced guiding- 
center Fokker-Planck collision operator ()47|) and the inclusion of effects associated with 
background sheared E x B velocities. 
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Appendix A. First-Order Guiding-Center Corrections 



The guiding-center transformation (x, p) — > (X^S, n,6) was originally derived by Lit- 



tlejohn ^20j by Lie-transform methods in the form of asymptotic expansions Z^^ 
eGf -l- ■ • ■, where the components of the first-order generating vector field are 



Gf 

G't 
G\ 



mc 

e 



2fi 



dpo 



mBitX: 



mvu 



B 



Po-R 



09 



b- Vb I - At (ai : Vb 
V1„B + Ja.:Vb+^ 



+ b- V X b 
(b • Vb • Po 



(A.l) 
(A.2) 
(A.3) 



and Gf = 0. Here, we use the rotating unit vectors (b, ±,p): p = b x ± = 1 cos6' — 
2 sin^, defined in terms of the fixed (local) unit vectors 1 X 2 = b, the vector field R = 
VI • 2 denotes Littlejohn's gyro-gauge vector field [20], which is used to define the gradient 
operator 



V* 



V + 



R + 



V X b b 



d_ 



(A.4) 



and the gyroangle- dependent dyadic matrices are defined as ai = — |(p± + ± p) = 032/06. 
Using expressions for and G^, the first-order correction to the gyroradius vector 



p^ IS given as 



Pi 



a Po X b - /3 Po - 7 b. 



(A.5) 
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where the coefficients a, /3, and 7 are defined as 



a 



/9 



7 



^ a, : Vb + ^ 

n ^ de 

( 1 "rfivl 
Or,- - V In 5 H ^ 



Vln5 



b-Vb , 



b-Vb + 



2fi 



b- V X b - ai : Vb 



9^ Vf] 



b- Vb 



^ib-Vln5 + 32 : Vb) . 



mfi2 V2 

Using Eq. ()A.5|) . the particle kinetic momentum ()41|) can now be expressed as 

b 



P + e 



X /X V5 + mwn'b- Vb 



^ ( da f II r ^ r 



90 



(A.6) 



where a, (3, and 7 are defined above. It is now a simple exercise using the formulas presented 
here to show that Eq. ()A.6|) satisfies property (ji^ . Lastly, the gyroangle averages used in 
Sec. IV. C [Eqs. (p^ - (jU7j) ] make use of the following expressions 



n dp. 



2 - -| b- V X b, 



P. • 



n dp. 



Vt dp. 



Pe • 



B 89 
_n_ dp, 

]Ib W 



4(^l-^b.Vxb), 



Pe 



2 pii b + (m V 



(A.7) 
(A.8) 
(A.9) 
(A.IO) 
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